Robust and efficient generator of almost maximal multipartite entanglement 
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Quantum chaotic maps can efficiently generate pseudo-random states carrying almost maximal 
multipartite entanglement, as characterized by the probability distribution of bipartite entanglement 
between all possible bipartitions of the system. We show that such multipartite entanglement is 
robust, in the sense that, when realistic noise is considered, distillable entanglement of bipartitions 
remains almost maximal up to a noise strength that drops only polynomially with the number of 
qubits. 

PACS numbers: 03.67.Mn, 03.67.Lx, 03.67.-a, 05.45.Mt 



Entanglement is not only the most intriguing feature 
of quantum mechanics, but also a key resource in quan- 
tum information science [l], Q • In particular, for quan- 
tum algorithms operating on pure states, multipartite 
(many-qubit) entanglement is a necessary condition to 
achieve an exponential speedup over classical computa- 
tion Q . The entanglement content of random pure quan- 
tum states is almost maximal 0, S 0| ! such states find 
applications in various quantum protocols, like super- 
dense coding of quantum states @,0|, remote state prepa- 
ration H, and the construction of efficient data- hiding 
schemes [9(. Moreover, it has been argued that random 
evolutions may be used to characterize the main aspects 
of noise sources affecting a quantum processor [loj . 

The preparation of a random state or, equivalently, the 
implementation of a random unitary operator requires 
a number of elementary one- and two-qubit gates ex- 
ponential in the number n q of qubits, thus becoming 
rapidly unfeasible when increasing n q . On the other 
hand, pseudo-random states approximating to the de- 
sired accuracy the entanglement properties of true ran- 
dom states may be gene rated efficiently, that is, polyno- 
Ll. Il2|. In particular, quantum chaotic 



^ ■ mially inn, ToL 11 



maps are efficient generators of multipartite entangle- 
ment among the qubits, close to that expected for ran- 
dom states 11, HI]. A related crucial question is whether 
the generated entanglement is robust when taking into 
account unavoidable noise sources affecting a quantum 
computer, that in general turn pure states into mixtures, 
with a corresponding loss of quantum coherence and en- 
tanglement content. In this paper we give a positive an- 
swer to this question. 

The number of measures needed to fully quantify 
multipartite entanglement grows exponentially with the 
number of qubits. Different measures capture various 
aspects of multipartite entanglement. Therefore, fol- 
lowing Ref. 14(, we characterize multipartite entangle- 



bipartitc entanglement between all possible bipartitions 
of the system. For pure states the bipartite entangle- 
ment is the von Neumann entropy of the reduced density 
matrix of one of the two subsystems: Eab{\^) ("01) = 
-Tr [p A log 2 pa] = S(p A ) , where p A = Tr B (\ip) (ip\), and 
A, B denote two subsystems made up of ua and Ub 
qubits {tia + ng = n q ). For sufficiently large systems 
[N = 2 rlq ^> 1), it is reasonable to consider only balanced 
bipartitions, i.e., with ua — wb, since the statistical 
weight of unbalanced ones becomes negligible Q . If the 
probability density has a large mean value (Eab) ~ n q 
(( • ) denotes the average over balanced bipartitions) and 
small relative standard deviation <tab/{Eab) <C 1, we 
can conclude that genuine multipartite entanglement is 
almost maximal (note that Eab is bounded within the 
interval [0, n g ]). This is the case for random states [14] . 

The model. The use of quantum chaos for efficient 
and robust generation of pseudo-random states carrying 
large multipartite entanglement is nicely illustrated by 
the example of the quantum sawtooth map 
map is described by the unitary operator U : 
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ment by means of a function rather than with a single 
measure: we look at the probability density function of 



where n = —id/ 88, [6,h] — i (we set h = 1) and the 
discrete time t measures the number of map iterations. 
In the following we will always study map ([T]) on the 
torus < 9 < 2ir, — tt < p < ir, where p = Tn. With 
an n g -qubit quantum computer we are able to simulate 
the quantum sawtooth map with N = 2 n « levels; as a 
consequence, 9 takes N equidistant values in the interval 
< 9 < 2tt, while n ranges from -N/2 to N/2 - 1 (thus 
setting T — 2-k/N). We are in the quantum chaos regime 
for map (p} when K = kT > or K < —4; in particular, 
in this work we focus on the case K = 1.5. 

There exists an efficient quantum algorithm for simu- 
lating the quantum sawtooth map [HI]. The crucial ob- 
servation is that the operator U in Eq. |T]) can be writ- 
ten as the product of two operators: Uu = e lk ^ e ^ 7 '^ I 2 
and Ut — e~' lTn / 2 , that are diagonal in the 9 and in 
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the n representation, respectively. Therefore, the most 
convenient way to classically simulate the map is based 
on the forward-backward fast Fourier transform between 
8 and n representations, and requires 0(N log 2 N) oper- 
ations per map iteration. On the other hand, quantum 
computation exploits its capacity of vastly parallelize the 
Fourier transform, thus requiring only 0((log 2 AO 2 ) one- 
and two-qubit gates to accomplish the same task [15] . In 
brief, the resources required by the quantum computer 
to simulate the sawtooth map are only logarithmic in 
the system size iV, thus admitting an exponential speed 
up, as compared to any known classical computation. 

Multipartite entanglement generation. We first com- 
pute (Eab) as a function of the number t of iterations of 
map ([TJ). Numerical data in Fig. [TJ exhibit a fast conver- 
gence, within a few kicks, of this quantity to the value 
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expected for a random state [J|. Precisely, {Eab) con- 
verges exponentially fast to (E™g d ), with the time scale 
for convergence oc n q (see inset of Fig. [TJ) . Therefore, the 
average entanglement content of a true random state is 
reached to a fixed accuracy within 0(n q ) map iterations, 
namely 0(n q ) quantum gates. We stress that in our case 
a deterministic map, instead of random one- and two- 
qubit gates as in Ref. [l2j, is implemented. Of course, 
since the overall Hilbert space is finite, the above expo- 
nential decay in a deterministic map is possible only up 
to a finite time and the maximal accuracy drops exponen- 
tially with the number of qubits. We also note that, due 
to the quantum chaos regime, properties of the generated 
pseudo-random state do not depend on initial conditions, 
whose characteristics may even be very different from it 
(e.g., in simulations of Fig. [TJ we start from completely 
disentangled states). 

As discussed above, multipartite entanglement should 
generally be described in terms of a function, rather than 
by a single number. We therefore show in Figf3]the prob- 
ability density function p(Eab) for the entanglement of 
all possible balanced bipartitions of the state |^t=3o)- 
This function is sharply peaked around (E™g d ), with a 
relative standard deviation tab/ (Eab) that drops expo- 
nentially with n q (see the inset of Fig. [2]) and is very small 
(~ 0.1) already at n q = 4. For this reason, we can con- 
clude that multipartite entanglement is large and that it 
is reasonable to use the first moment (Eab) of p(Eab) for 
its characterization. We have also calculated the corre- 
sponding probability densities for random states (dashed 
curves in Fig. [5]); their average values and variances are in 
agreement with the values obtained from states generated 
by the sawtooth map. The fact that for random states 
the distribution p(Eab) is peaked around a mean value 
close to the maximum achievable value £™Jf = n q /2 is 
a manifestation of the "concentration of measure" phe- 
nomenon in a multi-dimensional Hilbert space [B|, Q . 
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FIG. 1: (color online) Time evolution of the average bipartite 
entanglement of a quantum state, starting from a state of the 
computational basis (eigenstate of the momentum operator 
n), and recursively applying the quantum sawtooth map (fTJ) 
at K = 1.5 and, from bottom to top, n q = 4,6,8,10,12. 
Dashed lines show the theoretical values of Eq. ((2}. Inset: 
convergence of {EAB)(t) to the asymptotic value (Ea§ ). 
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FIG. 2: (color online) Probability density function of the bi- 
partite von Neumann entropy over all balanced bipartitions 
for the state \ipt), after 30 iterations of map ([TJ) at K — 1.5. 
Various histograms are for different numbers of qubits: from 



left to right n q 



10, 12; dashed curves show the corre- 



sponding probabilities for random states. Inset: relative stan- 
dard deviation (jab / (Eab) as a function of n q (full circles) 
and best exponential fit gab /{Eab) ~ e" ' 48 ™' (continuous 
line); data and best exponential fit gab /(Eab) ~ e~ n i /2 for 
random states are also shown (empty triangles, dashed line). 



Stability of multipartite entanglement. In order to as- 
sess the physical significance of the generated multipar- 
tite entanglement, it is crucial to study its stability when 
realistic noise is taken into account. Hereafter we model 
quantum noise by means of unitary noisy gates, that re- 
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suit from an imperfect control of the quantum computer 
hardware [l6| . We follow the noise model of Ref. 17 1. 
One-qubit gates can be seen as rotations of the Bloch 
sphere about some fixed axis; we assume that unitary 
errors slightly tilt the direction of this axis by a random 
amount. Two-qubit controlled-phase shift gates are diag- 
onal in the computational basis; we consider unitary per- 
turbations by adding random small extra phases on all 
the computational basis states. Hereafter we assume that 
each noise parameter Si is randomly and uniformly dis- 
tributed in the interval [—£, +e]; errors affecting different 
quantum gates are also supposed to be completely uncor- 
rected: every time we apply a noisy gate, noise param- 
eters randomly fluctuate in the (fixed) interval [— e, +e\. 

Starting from a given initial state IV'o); the quantum 
algorithm for simulating the sawtooth map in presence 
of unitary noise gives an output state \ip £l ^) that differs 
from the ideal output \ipt}- Here £j — (e\, e%, £ n4 ) 
stands for all the rid noise parameters Ei 1 that vary upon 
the specific noise configuration (nd is proportional to the 
number of gates). Since we do not have any a priori 
knowledge of the particular values taken by the param- 
eters £j, the expectation value of any observable A for 
our n g -qubit system will be given by Tr[/9 EjtJ 4], where the 
density matrix p £t is obtained after averaging over noise: 
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The integration over £/ is estimated numerically by sum- 
ming over Af random realizations of noise, with a statis- 
tical error vanishing in the limit Af — * oo. The mixed 
state p £ may also arise as a consequence of non-unitary 
noise; in this case Eq. ^ can also be seen as an unrav- 
eling of p £ into stochastically evolving pure states l^,), 
each evolution being known as a quantum trajectory [18j |. 

We now focus on the entanglement content of p £ j- 
Unfortunately, for a generic mixed state of n q qubits, 
a quantitative characterization of entanglement is not 
known, neither unambiguous (l9j . Anyway, it is possible 
to give numerically accessible lower and upper bounds 



for the bipartite distillable entanglement E^J(p £ ) 



max{S(p e , A ) - S(p £ ),0} < E ( ^(p £ ) < log 2 



= Tr 
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where p £ . A = Tr B (p £ ) and ||# 
notes the trace norm of the partial transpose of p £ 
respect to party B. 

In practice, we simulate the quantum algorithm for the 
quantum sawtooth map in the chaotic regime with noisy 
gates and evaluate the two bounds in Eq. Q for the dis- 
tillable entanglement of the mixed state p e j, obtained 
after averaging over Af noise realizations. A satisfac- 
tory convergence for the lower and the upper bound is 
obtained after Af ~ VN and Af ^ N noise realizations, 
respectively. In Fig. [3J upper panels, we plot the first mo- 
ment of the lower (E m ) and the upper (Em) bound for 




FIG. 3: (color online) Upper graphs: lower (E m ) (left panel) 
and upper bound (Em) (right panel) for the distillable entan- 
glement as a function of the noise strength at time t = 30. 
Various curves stand for different numbers of qubits: n q — 
4 (circles), 6 (squares), 8 (diamonds), 10 (triangles up), and 
12 (triangles down). Lower graphs: relative standard devia- 
tion of the probability density function for distillable entan- 
glement over all balanced bipartitions as a function of n q , 
for different noise strengths e. Dashed lines show a behavior 
a I (E) ~ e ~ n i/ 2 and are plotted as guidelines. 



the distillable entanglement as a function of the imperfec- 
tion strength. The various curves are for different num- 
bers n q of qubits; Af depends on n q and is large enough 
to obtain negligible statistical errors (smaller than the 
size of the symbols). In the lower panels of Fig. [3] we 
show the relative standard deviation of the probability 
density function (over all balanced bipartitions) for the 
distillable entanglement. Like for pure states, we no- 
tice an exponential drop with n q ; the distribution width 
slightly broadens when increasing imperfection strength 
e. We can therefore conclude that an average value of the 
bipartite distillable entanglement close to the ideal case 
e = implies that multipartite entanglement is stable. 

In order to quantify the robustness of multipartite en- 
tanglement with the system size, we define a perturbation 
strength threshold at which the distillable entangle- 
ment bounds drop by a given fraction, for instance to 
1/2, of their e = value, and analyze the behavior of 
£^ R ' as a function of the number of qubits. Numerical 
results are plotted in Fig. 0J both for lower and upper 
bounds we obtain a power-law scaling close to 
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It is possible to give a semi-analytical proof of the scal- 
ing (O for the lower bound measure, that is based on 
the quantum Fano inequality (20j . which relates the en- 
tropy S(p £ ) to the fidelity F = (ij>t\pe,t\$t)- S (P?) < 
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FIG. 4: (color online) Perturbation strength at which the 
bounds of multipartite entanglement halve (lower bound on 
the left panel, upper bound on the right panel), as a function 
of the number of qubits. Dashed lines are best power-law 



This result may become of practical relevance, since pro- 
totypes of quantum computers simulating these systems 
and, in particular, our specific model [22[ have been al- 
ready experimentally put on using a three-qubit NMR- 
based quantum processor 2^, 23[. The fact that distil- 
lable entanglement of balanced bipartitions remains al- 
most maximal, up to a noise strength which drops only 
polynomially with the number of qubits, supports the 
possibility that multipartite entanglement of a large num- 
ber of qubits might be used as a real physical resource in 
quantum information protocols. 
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h(F) + (1 - F) log 2 (iV 2 -l), where h{x) = -xlog 2 (x) - 
(1 — x) log 2 (l — x) is the binary Shannon entropy. Since 
F ~ e -^ 2n s' [IE ED], with 7 - 0.28 and n g = 3n 2 + n q 
being the number of gates required for each map step, we 
obtain, for e 2 n„t <C 1, 



S(p s ) < ie 2 n g t 



log 2 (7£ 2 ?i 5 i) + 2n q 



1 
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For sufficiently large systems the second term dominates 
(for n q — 12 qubits, t — 30 and e~5x 10 -3 the other 
terms are suppressed by a factor ~ 1/10) and, to a first 
approximation, we can only retain it. On the other hand, 
an estimate of the reduced entropy S(pa) is given by the 
bipartite entropy (J2J> of a pure random state There- 
fore, from Eq. ((4]) we obtain the following expression for 
the lower bound of the distillablc entanglement: 
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From the 
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threshold definition E A D ^{p e ( R )) 
^S(pa) we get the scaling (J3J) , that 
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is valid when n q 3> 1: e 



l/y/247 n 2 t. Notice that, 

for small systems as the ones that can be numerically 
simulated (see data in Fig. , the first term of Eq. ^ 
may introduce remarkable logarithmic deviations from 
the asymptotic power-law behavior. At any rate, the 
scaling derived from Eq. is in good agreement with 
our numerical data, and also reproduces the prefactor in 
front of the power-law decay ([5]) up to a factor of two. 

Conclusions. We have shown that quantum chaotic 
maps provide a convenient tool to efficiently generate 
in a robust way the large amount of multipartite entan- 
glement close to that expected for truly random states. 
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